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ABSTRACT 

A  theoretical  model  of  expendable  bathythermograph  (XBT)  fall  rate  is  reviewed,  and  a  new  form  of  fall- 
rate  equation  is  proposed  to  include  near-surface  transient  effects.  Comparisons  are  made  of  T-7  XBT  and  CTD 
(conductivity,  temperature,  and  depth )  depths  of  thermohaline  features  off  Barbados.  Fall-rate  equation  coef¬ 
ficients  are  derived  and  compared  with  the  manufacturer-supplied  coefficients.  As  other  investigators  have 
found,  the  Sippican  equation  consistently  underestimates  probe  depth  by  as  much  as  35  m  at  760  m.  Analysis 
yields  a  new  equation,  Z  =  6.798/  -  0.002383/2  -  4.01,  for  depths  greater  than  about  10  m.  Considerable 
probe-to-probe  variability  is  noted  and  is  found  to  be  primarily  the  result  of  differences  in  the  linear  term  or 
terminal  velocity  of  the  probes;  variation  in  effective  drag  resulting  from  probe  irregularities  is  the  likely  cause. 
Recommendations  for  additional  work  are  made. 


1.  Introduction 

The  expendable  bathythermograph  (XBT)  was  in¬ 
troduced  by  the  Sippican  Corporation  in  1965.  Since 
then  it  has  become  a  widely  used  method  of  determin¬ 
ing  temperature  profiles  in  the  ocean,  and  is  often  de¬ 
pended  on  as  a  primary  observational  tool.  There  has 
been  recurrent  interest  in  improving  the  accuracy  of 
XBT  measurements,  particularly  the  determination  of 
probe  depth.  Since  XBT  depth  is  not  measured  directly 
but  inferred  from  elapsed  time  and  a  fall-rate  model, 
a  number  of  factors  affect  the  accuracy  of  the  result. 
That  these  factors  are  variable  has  been  documented 
by  several  investigators  whose  results  have  been  sum¬ 
marized  by  Singer  ( 1990)  in  his  Table  3.  In  the  present 
study  we  examine  the  fall-rate  model  and  the  variability 
of  factors  affecting  it  with  a  suite  of  comparison  data. 

The  depths  of  most  expendable  probes  [i.e.,  XBTs, 
expendable  sound  velocimeters  (XSVs),  expendable 
current  profilers  (XCPs)]  are  usually  calculated  using 
an  equation  of  the  form 

Z  =  bt  —  at2,  ( 1 ) 

where  t  is  time  elapsed  from  probe  entry  into  the  water. 
Since  about  1974  it  has  been  recognized  that  this  equa¬ 
tion  has  some  inadequacies,  particularly  in  the  values 
of  the  coefficients.  Singer  ( 1990)  provides  a  concise 
summary  of  past  evaluations  of  Eq.  ( 1 )  for  the  T-7 
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XBT.  In  all  cases  it  was  found  that  the  manufacturer- 
supplied  coefficients  ( a  =  0.00216  m  s-2,  b  -  6.472 
m  s-1 )  are  systematically  too  small. 

In  the  present  study  we  examine  T-7  XBT  and  con¬ 
current  CTD  (current,  temperature,  and  depth)  tem¬ 
perature  profiles  acquired  in  a  region  near  Barbados, 
where  the  thermohaline  structure  of  the  upper  water 
column  is  ideal  for  such  comparisons.  The  purpose  of 
the  study  is  to  further  examine  the  validity  of  Eq.  ( 1 ) 
for  determining  XBT  depth.  In  section  2  we  review 
the  theoretical  basis  for  Eq.  ( 1 )  presented  by  Green 
( 1984),  hereafter  G84,  and  suggest  a  modification  to 
the  fall-rate  equation  (FRE).  Section  3  is  a  brief  syn¬ 
opsis  of  previous  studies  of  the  XBT  fall-rate  problem. 
In  section  4  we  describe  field  measurements  and  present 
an  analysis  of  the  resulting  XBT  and  CTD  data  and 
discuss  implications  for  the  FRE.  Finally,  in  section  5, 
we  make  some  recommendations  for  additional  work 
that  will  be  required  if  real  improvements  to  expend¬ 
able-probe  depth  accuracies  are  to  be  made. 

2.  XBT  fall  rates:  A  simple  model 

While  processes  involved  in  the  descent  of  an  ex¬ 
pendable  probe  may  be  quite  complex  in  detail,  a  rel¬ 
atively  simple  analytical  model  can  describe  most  of 
the  observational  results.  The  model  is  that  described 
by  G84,  and  is  presented  here  in  a  somewhat  simplified 
way.  A  balance  is  assumed  involving  the  vertical  ac¬ 
celeration,  the  net  buoyant  force,  and  the  hydrody¬ 
namic  drag,  the  last  being  proportional  to  the  square 
of  the  speed, 

d2Z  g(m  -  mw)  pHCDAW2 

dt 2  m  2m  ’  {  ’ 
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where  Z  is  vertical  position  ( positive  downward )  and 

pH  =  water  density 
Cp  =  drag  coefficient 
A  =  effective  cross-sectional  area  of  probe 
W  =  vertical  speed 
m  =  mass  of  probe 
m»  =  mass  of  displaced  water 
g  =  acceleration  of  gravity 
t  =  time. 

The  parameters  in  Eq.  ( 2 )  can  be  expected  to  vary 
from  region  to  region  (pM  ,  m„,  CD)  and  from  probe  to 
probe  (m,CD,A).  Furthermore,  it  is  clear  that  m  varies 
with  depth  (as  expendable  probe  wire  unreels),  and  it 
is  likely  that  CD ,  m„  ,  and  p„  will  vary  with  depth,  as 
well.  Equation  ( 1 )  implicitly  assumes  that  the  left  side 
of  Eq.  (2)  is  negligible;  here  fall  speed  is  assumed  to 
be  virtually  equal  to  the  terminal  velocity,  although 
the  latter  may  change  slowly  with  depth.  This  assump¬ 
tion  is  reasonable  for  depths  in  excess  of  about  10  m 
(G84),  and  we  define  this  balance  as  asymptotic  case 
1.  In  this  case, 

(J) 

where  m'  =  m  -  mw.  The  parameters  p„,  CD,  and  m' 
can  vary  with  depth.  The  last  of  these  varies  as  a  result 
of  loss  of  XBT  wire  as  the  probe  falls.  The  variation 
of  W2  due  to  the  change  in  water  density  with  depth 
is  at  most  about  10%  of  that  due  to  wire  payout,  and 
will  be  neglected  here.  The  drag  coefficient  also  can 
vary  ;  the  character  and  magnitude  of  this  change  is  as 
yet  a  matter  of  controversy  (G84). 

Wire  payout  is  effectively  modeled  as  a  linear  func¬ 
tion  of  depth: 

m'  =  mo(  1  -  tZ).  (4) 

Based  on  G84  and  discussions  with  the  manufacturer 
we  estimate  (  =  1.29  X  '0-4  nT1 .  Following  G84  we 
also  model  the  variation  in  CD  as  a  linear  function  of 
depth: 

CD  =  CD0(\  +SZ).  (5) 

Then,  taking  the  linear  terms  of  an  expansion. 


-2(1  -ijZ)' 


=  WTt  +  C. 


If  we  take  Z(0)  =  0,  then  C  =  -2/tj,  and 

z. (9) 

4 

which  is  the  result  obtained  by  G84  and  is  the  same 
form  as  Eq.  ( 1 ).  The  quadratic  term  depends  on  the 
basic  (initial)  terminal  velocity  and  on  the  parameter 
t)  (which  depends  primarily  on  the  probe’s  mass  loss 
due  to  wire  payout).  Using  the  manufacturer-supplied 
coefficients  we  can  estimate  ij: 

r^CL=  0.00216 
4 

WT  =  6.472, 

so  t ?  =  2.06  X  10~4.  Then  from  Eq.  (6) 

5  =  t)  -  t  =  (2.06  -  1.29)  X  10'4  =  7.7  X  10“5. 

For  Z  =  10m(/«=  1.5  s),  the  contribution  of  the 
quadratic  term  is  only  about  0.5  cm.  Hence,  for  this 
shallow  depth  interval  we  examine  asymptotic  case  2, 
where  we  ignore  variations  in  parameters,  but  do  in¬ 
clude  the  acceleration  term.  We  rewrite  Eq.  (2)  as 


where 


PwCDA 


\aj  \pwCDA ) 


(as  before).  Rearranging: 


W  f  -  W‘ 


—  adt. 


( 1  —  fZ  -  5Z)  =  W  j-(  1  -tjZ), 

PwCpoA 

(6) 

where  rj  =  <  +  8.  This  can  be  directly  integrated  to  yield 
depth  of  the  probe  as  a  function  of  time, 

+  (7) 

where  W\=  2gm’n/ pwCnoA .  Upon  evaluation  of  the 
integral , 


The  integral  of  this  is 


Note  the  two  regimes  that  depend  on  the  ratio  W/WT 
(there  is  a  singularity  for  W  =  fFr). 

Applying  the  initial  condition  at  t  =  0  where  W 
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=  W0  we  evaluate  C.  After  some  rearranging  we  then  Eq.  ( 1 )  is  negligible.  Hence,  it  is  reasonable  to  propose 
have  a  combined  FRE  as 


W=  WT 


coth(al*V/  +  C,);  -^>1 

W  f 

tmh(aWTt  +  C2);  -^<1 

wT 


,  (13) 


where  Cj  =  co\\i~x  (W0IWT)  and  C2  =  tanh“‘(Ifo/ 
WT).  Equation  (13)  is  then  integrated  to  find  Z(t). 
After  applying  the  condition  Z(0)  =  0.  we  obtain 


r 

'sinhfaW'Vf  +  C\ )" 

•  i^>i] 

’  WT 

1 

In 

sinhCj 

- 

a 

In 

'cosh(aWV  +  C2) 

W0 

coshC2 

'  WT  < 

(14) 


Table  1  lists  W,  Z,  and  A Z  =  Z  -  WTt  for  three  cases 
of  WQ.  The  largest  absolute  values  of  A Z  are  seen  for 
W0  =  0  and  suggest  an  offset  of  about  4  m  to  be  applied 
deeper  than  about  10  m  when  using  Eq.  ( 1 )  to  calculate 
depth.  If  H  o  is  close  or  equal  to  H V  this  transient  effect 
is  minimal  and  A Z  will  be  close  to  zero.  The  T-7  XBTs 
are  typically  dropped  from  about  2  m  above  the  sea 
surface,  which  should  ( assuming  a  clean,  vertical  entry) 
result  in  an  entry  velocity  close  to  the  nominal  value 
of  Wr,  about  6.5  m  s_1.  Other  factors  not  modeled, 
such  as  spinup  of  the  probe,  non  vertical  entry,  and  air 
entrapped  in  the  wire  spool  may  also  slow  the  probe 
in  the  first  few  meters. 

It  can  be  shown  that  for  t  sufficiently  large  (about 
1.5  s  for  nominal  parameters)  such  that  Z  is  of  the 
order  of  10  m  or  larger,  Eq.  ( 14)  approaches  the  form 

Z  =  WTt  +  D, 


where  D  (the  offset )  depends  on  the  initial  speed,  W0. 
Here  D  is  negative  for  W0  <  WT  and  positive  for  the 
reverse.  For  W0  =  WT,  D  =  0.  We  have  shown  that 
for  t  of  the  order  of  1 .5  (or  less)  the  quadratic  term  in 


Z  =  Zdt, 


(15) 


where  Z(  is  just  Eq.  ( 14).  A  relationship  between  W0 
and  D  can  be  derived  from  Eq.  (14),  valid  for  all  pos¬ 
itive  values  of  the  ratio  W0/WT\ 

TTT  =  2eaD  -  1,  (16) 

or 

„  1,  (lTo/BV)+l 

D  =  -  In - - - . 

a  2 


Hence,  an  apparent  W0,  based  on  an  empirically  de¬ 
termined  value  of  D,  can  be  used  in  Eq.  ( 15)  to  de¬ 
termine  Z.  The  minimum  value  for  D  (for  a  given  a) 
is  just 


T^min 


ln(0.5)  -0.6931 


a 


(17) 


which  occurs  for  W0  =  0. 


3.  Previous  XBT  fall-rate  studies 

Attempts  to  improve  the  XBT  FRE  are  not  new. 
McDowell  ( 1977),  Fedorov  et  al.  ( 1978),  and  Flierl 
and  Robinson  ( 1977)  discuss  systematic  offsets  of  iso¬ 
therm  depths  derived  from  XBT  and  CTD  data.  The 
errors  they  found  were  generally  of  the  same  form  and 
magnitude.  This  suggested  a  possible  error  in  the  man¬ 
ufacturer’s  FRE  coefficients.  Seaver  and  Kuleshov 
( 1 982 )  proposed  a  somewhat  complex  scheme  to  cor¬ 
rect  XBT  depths.  They  considered  the  effects  of  tem¬ 
perature  and  salinity  on  drag.  Heinmiller  et  al.  ( 1983) 
also  observed  this  systematic  error  and  proposed  a 
scheme  to  correct  T-7  XBT  depths.  This  was  a  piece- 
wise  linear  correction  to  the  recorded  depth,  with  two 
segments  meeting  at  325  m;  the  error  profile  ( i.e.,  ZXbt 


Table  1.  Depth  [From  Eq.  (14)].  fall  rate  [from  Eq.  (13)],  and  AZ  =  Z  -  WTt  for  three  values  of  H  o,  a  =  0.1835  m*', 
0  =  7.675  m  s-2,  and  WT  =  6.467  m  s"1  (parameter  values  adapted  from  G84). 


Time 

(s) 

Ho/HV  =  0 

W<JWT  =  0.5 

H'o/H V  -  2 

Z 

W 

A  Z 

Z 

W 

AZ 

Z 

W 

AZ 

0.00 

0.00 

0.000 

0.00 

0.00 

3.234 

0.00 

0.00 

12.936 

0.00 

0.40 

0.59 

2.859 

-1.99 

1.68 

4.990 

-0.91 

4.04 

8.384 

1.46 

0.80 

2.16 

4.783 

-3.02 

3.87 

5.853 

-1.30 

7.11 

7.148 

1.93 

1.20 

4.29 

5.759 

-3.47 

6.30 

6.223 

-1.46 

9.87 

6.723 

2.10 

1.60 

6.69 

6.184 

-3.66 

8.82 

6.372 

-1.53 

12.52 

6.566 

2.17 

2.00 

9.21 

6.357 

-3.73 

11.38 

6.431 

-1.55 

15.13 

6.506 

2.19 

2.40 

11.76 

6.425 

-3.76 

13.96 

6.454 

-1.56 

17.73 

6.483 

2.20 

2.80 

14.34 

6.451 

-3.77 

16.54 

6.462 

-1.57 

20.32 

6.474 

2.21 

3.20 

16.92 

6.462 

-3.78 

19.13 

6.466 

-1.57 

22.91 

6.470 

2.21 

3.60 

19.51 

6.466 

-3.78 

21.72 

6.467 

-1.57 

25.49 

6.469 

2.21 
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-  ZctD)  appeared  to  be  roughly  linear  in  the  two  re¬ 
gimes.  They  suggested  that  some  of  the  error  may  have 
been  due  to  the  analog  recordings  used  at  the  time. 
Hanawa  and  Yoritaka  (1987)  discussed  FRE  coeffi¬ 
cient  modification  based  on  data  from  Japanese-made 
XBTs  and  a  digital  acquisition  system.  They  used  a 
regression  analysis  that  resulted  in  coefficients  larger 
than  those  of  the  manufacturer.  The  error  profile  they 
found  was  nearly  linear,  unlike  those  found  previously. 

Wright  and  Szabados  ( 1 989 )  examined  the  FRE  for 
the  T-7  (among  others)  with  data  acquired  in  the  same 
region  as  the  present  study.  They  also  found  that  the 
manufacturer’s  FRE  underestimates  depths  of  thermal 
features  observed  with  concurrent  CTD  casts.  They 
used  a  regression  scheme  that  adjusted  the  apparent 
sampling  time  interval  of  the  XBTs  to  minimize  the 
error  profile.  The  procedure  led  to  correction  factors 
for  the  FRE  coefficients.  Their  results  differed  from 
those  of  Hanawa  and  Yoritaka  ( 1987)  by  about  1 1  m 
at  a  depth  of  750  m.  Since  these  datasets  were  acquired 
in  different  oceans  with  probes  manufactured  in  dif¬ 
ferent  countries,  some  differences  should  be  expected. 

Singer  (1990)  presented  a  summary  of  previous 
studies  of  the  T-7  FRE  and  compared  them  to  results 
using  data  from  the  Gulf  of  Mexico.  His  FRE  coeffi¬ 
cients  fall  between  those  found  by  Hanawa  and  Yori¬ 
taka  ( 1987)  and  Wright  and  Szabados  ( 1989).  Singer 
also  found  an  offset  term  (about  4  m)  for  his  quadratic 
FRE;  he  stated  that  this  offset  was  “.  .  .  an  artifact  of 
the  linear  fit  .  .  .”  which,  unlike  other  investigators’ 
equations,  had  not  been  found  to  fit  to  zero  depth  at 
zero  time.  Indeed,  we  infer  that  the  fall  rate  is  not  linear 
close  to  the  surface,  as  is  also  suggested  by  Eq.  ( 15). 

4.  Observations 

In  this  section  we  briefly  describe  (a)  a  test-and- 
evaluation  experiment  involving  CTD  and  XBT  mea¬ 
surements.  We  then  describe  (b)  comparisons  and 
analyses  of  these  data  and  the  determination  of  new 
FRE  coefficients  using  selected  features  in  the  temper¬ 
ature  profile,  and  we  subsequently  discuss  (c)  the  vari¬ 
ability  of  individually  determined  (for  each  XBT 
probe)  FRE  coefficients.  Finally  (d)  we  examine  the 
implications  of  the  new  FRE  for  the  entire  XBT  range. 

a.  Data 

In  May  1990  we  conducted  CTD-expendable-probe 
comparisons  near  16°N,  56°W,  about  300  nm  north¬ 
east  of  Barbados.  The  thermohaline  structure  of  the 
upper  100  m  in  this  region  consists  of  large,  well-de¬ 
fined,  and  persistent  “staircases”  (alternating  isother¬ 
mal,  isohaline  layers  and  high-gradient  “sheets”). 
These  characteristics  are  ideal  for  comparisons  of  hy¬ 
drographic  profiling  systems.  These  structures  are  par¬ 
ticularly  useful  in  the  present  case  because  depth  un¬ 
certainties  can  be  effectively  decoupled  from;  for  ex¬ 


ample,  temperature  uncertainties.  In  this  paper  we 
describe  a  subset  of  the  data  acquired  during  the  overall 
experiment. 

Over  a  period  of  about  2  h,  36  sets  of  4  simultaneous 
T-7  XBTs  were  dropped  using  four  Sippican  Mark-9 
Launcher  Acquisition  Systems  (LAS’s).  Concurrently, 
five  CTD  casts  were  made  to  about  1000  m  with  a  Neil 
Brown,  Mark  III  CTD  system.  During  this  period,  the 
ship  was  dead  in  the  water  and  conditions  were  good. 

The  CTD  was  calibrated  at  the  Naval  Oceanographic 
Office  before  the  cruise.  Temperature  error  was  found 
to  be  less  than  0.003°C,  and  pressure  uncertainty  was 
less  than  4  db  between  0  and  1 000  db.  Additionally,  a 
short  time  series  ( about  1 0  min )  of  CTD  data  collected 
with  the  probe  at  a  measured  average  depth  of  4  m 
resulted  in  an  average  pressure  of  2.5  db,  indicating 
that  the  CTD  pressure  was  too  low  by  1.5  db.  CTD 
pressures  were  corrected  by  this  amount  prior  to  further 
analysis,  and  converted  to  depth  by  integrating  the 
CTD- measured  specific  volume  of  seawater  using  a 
constant  value  of  the  acceleration  of  gravity  of  9.8 
m  s  ’.  Of  the  144  probes  (T-7’s)  dropped  between 
about  0545  and  0800  LJTC  only  128  were  recorded. 
Data  from  16  probes  were  lost  as  a  result  of  acquisition 
system  malfunction;  most  of  these  ( 14)  losses  occurred 
on  LAS  3.  Of  the  recorded  data,  10  probes  (8%)  pro¬ 
duced  bad  temperature  data  or  ended  early  (broken 
wire).  A  typical  set  of  simultaneous  XBT  temperature 
profiles  is  shown  in  Fig.  I.  Superimposed  (dashed 
curve)  is  a  concurrent  (time-interpolated)  CTD  tem¬ 
perature  profile.  Since  CTD  and  XBT  temperatures 
were  sampled  at  different  times  at  each  depth  due  to 
the  sampling  scheme  and  widely  different  descent  rates 
and  since  the  thermal  field  was  changing  on  short  pe¬ 
riods  due  probably  to  internal  wave  straining,  time  in¬ 
terpolation  of  the  CTD  data  was  necessary.  In  partic¬ 
ular,  depth-interpolated  CTD  profiles  were  interpolated 
in  time  to  coincide  with  XBT  sample  times  (the  time- 
interpolation  procedure  is  further  explained  next). 

Two  things  are  immediately  evident  in  Fig.  1 .  The 
CTD  data  show  the  step  features  significantly  deeper 
than  do  the  XBT  data,  consistent  with  comparisons 
made  previously  ( Hanawa  and  Yoritaka  1987;  Wright 
and  Szabados  1989;  and  Singer  1990).  The  simulta¬ 
neously  dropped  XBTs  show  differences  in  feature 
depths  among  themselves.  The  latter  observation  holds 
to  a  greater  or  lesser  extent  (than  the  case  shown  in 
Fig.  I )  for  all  simultaneous  sets  of  XBTs. 

b.  Regressions  based  on  selected  feature  depths 

To  better  quantify  the  depth  discrepancies,  expanded 
plots  of  segments  of  the  profiles  containing  identifiable 
features  were  examined.  For  the  example  (Figs.  2a,b) 
depth  discrepancies  between  simultaneous  XBT  pro¬ 
files  exhibit  both  overall  offsets  and  depth  dependen¬ 
cies.  Also  evident  are  temperature  offsets  in  the  iso¬ 
thermal  layers  (up  to  about  0.1  °C  in  the  example). 
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Temperature(deg  C) 

Fig.  1.  Typical  set  of  simultaneous  XBT  profiles  (solid  curves) 
and  time-interpolated  CTD  profile  (dashed)  of  temperature. 


Most  such  temperature  offsets  in  these  data  were  within 
the  manufacturer’s  tolerance  of  0.2°C,  but  they  made 
automated  analysis  algorithms  difficult  to  formulate 
(some  of  the  temperature  errors  were  functions  of 
depth).  Hence,  we  decided  to  identify  manually  a  set 
of  features  that  were  common  to  all  XBTs  and  CTDs. 
Figures  3a,b  show  the  16  features  selected  on  one  of 
the  example  XBTs.  The  depths  of  these  points  were 
manually  digitized  with  a  resulting  estimate  of  ±1  m 
digitization  error.  Rather  than  digitizing  the  time-in¬ 
terpolated  CTD  profiles  (such  as  that  shown  in  Figs. 
1  and  2)  the  five  original  CTD  temperature  profiles 
were  digitized  and  the  feature  depths  subsequently  in¬ 
terpolated  to  respective  XBT  times.  The  time-inter¬ 
polation  scheme  assumed  an  average  CTD  descent  rate 
of  1.02  ±  0.04  m  s_1  (based  on  actual  CTD  profiles) 
yielding  an  uncertainty  in  the  time  of  interpolated  Zero 
of  about  1  min  (the  same  scheme  was  used  to  time 
interpolate  the  full  CTD  profiles  discussed  previously). 
XBT  data  were  recorded  as  depth  and  temperature  so 
the  time  of  each  XBT  sample  was  calculated  by  in¬ 
verting  the  manufacturer’s  FRE: 

Ixbt  =  1498 [  1  -  0.5(1  -  2.0627  X  10-4ZXBt),/2]. 

(18) 

where  Z  is  the  depth  in  meters  and  t  is  in  seconds. 
Time  interpolation  provided  a  CTD  depth  for  each 


1femperature(deg  C) 


Fig.  2.  Typical  set  of  simultaneous  XBT  profiles  (solid  curves) 
and  time-interpolated  CTD  profile  (dashed)  of  temperature:  (a)  upper 
140  m,  (b)  400-700  m. 
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Tempera turefdeg  C) 


Fig.  3.  Typical  XBT  profile  showing  selected  features  (dots)  used 
in  the  analysis:  (a)  upper  140  m,  (b)  400-700  m. 


XBT  group,  /,  and  each  feature,,/.  Similarly,  for  each 
feature  and  XBT  group  we  have  a  set  of  up  to  four 
XBTs,  k.  Hence,  ZctdUJ)  is  a  function  of  (xbtO,  V, 
k)  for  /  =  1,  36;  j  =  1,  16;  k  =  1,  km!l%  4,  where  we 
take  Zctd  to  be  the  “correct”  feature  depth,  and  XBT 
depths  with  their  errors  are  expressed  as  (Xbt  through 
Eq.(18). 

Depth  errors  (Zero  -  ZXbt  )  for  all  features  are  plot¬ 
ted  as  functions  of  Zctd  in  Fig.  4.  An  offset  of  ap¬ 
proximately  4  m  near  the  surface  and  the  positive  bias 
at  depth  reported  by  other  investigators  are  apparent. 
Additionally,  the  large  scatter  suggested  by  the  profile 
plots  is  greater  at  depth  than  near  the  surface.  Of  these 
differences,  67%  exceed  the  manufacturer’s  cited  tol¬ 
erance  of  2%  of  depth  or  5  m. 

The  manufacturer’s  FRE  [Eq.  ( 1 )]  is  quadratic  in 
time  but  contains  no  constant  term;  it  is  of  the  same 
form  as  Eq.  (9)  for  asymptotic  case  1.  To  account  for 
possible  offsets  resulting  from  transient  effects  as  in 
asymptotic  case  2  or  some  other  systematic  data  offset, 
we  included  an  unspecified  constant  term  and  sought 
all  three  coefficients  by  least-squares  analysis.  That  is, 
we  minimized 

E2  =  (Z-  Z)2,  (19) 

where 

Z=Zcn>(i,j )  (20) 

Depth  Error  (m) 

-25  -20  -15  -10  -5  0  5  10  15  20  25  30  35  40 


Fig.  4.  Depth  error  (Zctd  -  ZXbt)  for  selected  features.  XBT 
depths  calculated  with  the  Sippican  FRE.  Different  symbols  denote 
different  LAS's. 
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Table  2.  Fall-rate  equation  coefficients  from  regression  analysis  with  uncertainty  error  e  and  overall  rms  error  of  fit  E. 


a  (m  s  J) 

e{%) 

b( ms  ') 

e(%) 

c(m) 

e(%) 

E(  m) 

Case  1  (fully  quadratic) 

0.00409 

5.7 

6.962 

4.1 

5.18 

Case  2  (modified) 

0.00238 

— 

6.798 

5.25 

and 

Z  =  k)  +  btXBT(i,j,  k)  +  c.  (21) 


The  overbar  denotes  an  average  over  a  set  of  Z,  /  pairs. 

In  the  first  case  we  averaged  over  all  1 1 8  good  profiles 
and  all  16  features  in  each  profile  (a  1888-member 
average).  Coefficients  resulting  from  a  straightforward 
regression,  with  associated  errors,  appear  in  Table  2. 
From  Eq.  (9)  in  section  1  we  derived  from  these  a 
value  for  r),  the  linear  component  of  depth  dependence 
of  fall-rate  parameters: 


4,102  X  10~3  X  4 
(6.95)2 


=  3.397  X  10~4  nT1. 


(22) 


This  is  about  1 .6  times  the  value  implied  by  the  original 
FRE.  Since  t?  =  t  +  5  and  e  =  1.29  X  10“4  m_1  [Eq. 
(4)],  5  =  2.1  X  10  4  m_l,  an  unreasonably  large  change 
of  the  drag  coefficient  (G84)  is  implied.  Furthermore, 
the  quadratic  and  linear  coefficients  are  larger  than 
those  found  by  Wright  and  Szabados  ( 1989)  for  the 
same  region. 

An  alternative  approach  constrains  the  relationship 
between  a  and  b  by  specifying  a  plausible  value  of  r\. 
Since  we  do  not  really  know  6  exactly,  we  use  the  value 
of  7]  inferred  from  the  original  FRE,  2.06  X  10  ~4  m_l . 
Here,  a  modified  regression  was  performed  using  a 
specified,  initial  value  for  a.  The  coefficients  b  and  c 
were  found,  and  a  new  value  of  a  was  calculated  using 
the  fixed  r;  and  Eq.  (22).  This  process  was  repeated 
until  the  coefficients  converged  to  constant  values 
(typically  abut  four  iterations).  These  also  appear  in 
Table  2  (no  uncertainty  is  presented  for  a  here,  since 
it  is  specified  and  dependent  on  b).  The  rms  error  in 
this  case  is  only  slightly  greater  than  the  error  from  the 
previous  approach.  Furthermore,  the  quadratic  and 
linear  coefficients  are  virtually  identical  to  those  found 
by  Wright  and  Szabados  ( 1989).  Indeed,  in  their  anal¬ 
ysis,  while  apparently  different  from  that  presented 
here,  the  quadratic  and  linear  terms  were  linked;  they 
used  only  one  free  parameter,  a  sampling  interval  cor¬ 
rection  (they  did  not  allow  a  constant  offset).  Hence, 
they  effectively  constrained  their  analysis  to  the  original 
value  of  t),  as  has  been  done  here. 

The  Li  Terences  between  depths  calculated  with  the 
Sippican  FRE  and  those  calculated  with  the  two  FREs 
represented  in  Table  2,  as  well  as  FREs  proposed  by 
other  investigators,  are  plotted  as  functions  of  depth 
in  Fig.  5.  The  curve  of  our  second  FRE  (Table  2)  is 


parallel  to  that  of  Wright  and  Szabados  ( 1989),  dif¬ 
fering  only  by  the  offset  term.  Singer’s  (1990)  equations 
[(3a)  and  (3b)  of  his  Table  3]  appear  to  yield  quite 
similar  results;  his  equations  differ  by  2  m  at  most  ( at 
760  m),  which  is  insignificant  in  the  light  of  the  un¬ 
certainties  we  found.  Singer’s  Eq.  (3a),  the  equations 
proposed  by  Hanawa  and  Yoritaka  ( 1987 )  and  Wright 
and  Szabados,  and  our  second  equation  exhibit  nearly 
linear  difference  curves,  differing  only  slightly  in  slope 
and  offset.  Our  first  equation  is  clearly  the  anomaly 
because  of  the  relatively  large  quadratic  coefficient.  Our 
two  equations  agree  at  65  and  56 1  m.  There  is  a  relative 
maximum  discrepancy  of  about  2.5  m  near  a  depth  of 
320  m  and  a  maximum  discrepancy  at  the  bottom  of 
the  profile  (760  m)  of  about  6  m  (the  disagreement 
between  the  derived  equations  and  the  Sippican  equa¬ 
tion  exceeds  25  m  at  760  m) .  With  the  rms  error  of  fit 
for  both  derived  equations  greater  than  5  m,  there  is 
no  clear  statistical  basis  for  selecting  one  or  the  other 
of  these  sets  of  coefficients  as  the  more  correct.  Since 
the  second  case  constrains  the  quadratic  and  linear 
coefficients  to  be  consistent  with  a  reasonable  physical 
model,  we  selected  this  case  for  all  subsequent  discus¬ 
sion.  Future  investigations  should  include  efforts  to  re¬ 
fine  estimates  of  the  parameter  77. 

Singer  ( 1990)  summarized  XBT  fall-rate  studies  by 
other  investigators.  He  found  a  similar  negative  offset 
(about  4  m)  to  that  found  here  (the  others  apparently 
did  not )  and  suggested  that  the  probes  may  have  fallen 
more  slowly  between  the  surface  and  100  m  than  they 
did  deeper  than  this  level.  An  apparent  lower  descent 
rate  near  the  surface,  hence  a  negative  offset  term  in 
the  regression,  would  result  if  the  probe  were  to  begin 
its  descent  at  less  than  its  terminal  velocity.  We  at¬ 
tempted  to  reconcile  the  observed  offset  term  with  the 
model  described  in  section  2.  A  plausible  value  for  a 
[in  Eq.  (17)]  is  0.185  m_l  yielding  Dm in  =  -3.7  m. 
Hence,  if  the  model  implicit  in  Eq.  (15)  is  relevant, 
W0  would  have  to  be  zero.  Since  the  probes  were 
dropped  from  a  height  of  1.5-2  m  above  the  surface. 
H’0  was  probably  clo«p  to  W >,  suggesting  that  the  model 
may  be  inadequate  or  that  there  may  have  been  initial 
decelerations  when  the  probes  impacted  the  water. 

Using  the  coefficients  resulting  from  the  modified 
regression  analysis.  XBT  depths  were  recalculated. 
These  depth  errors  for  the  features  are  plotted  in  Fig. 
6.  Now,  87%  of  the  points  fall  within  the  cited  tolerance 
(indicated  by  the  lines  on  the  plot).  A  similar  calcu¬ 
lation  using  the  coefficients  from  the  fully  quadratic 
case  produced  a  very  similar  distribution  ( not  shown ). 
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Depth  Difference  (m) 
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Fig.  5.  Differences  between  depth  calculated  with  the  Sippican  FRE  and  other  FREs  discussed 
in  the  text.  Specifically  ZSippK.,n  -  Z,,  where  i  =  HT-1  (first  equation  of  Table  2),  HT-2  (second 
equation  of  Table  2),  HY  (Hanawa  and  Yoritaka  1989),  WS  (Wright  and  Szabados  1987), 
Sg-a  [Singer  1990,  Eq.  (3a)],  Sg-b  [Singer,  Eq.  (3b)]. 


Depth  Error  (m) 
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Fig.  6.  Depth  error  (Z<TD  -  ZXBT)  for  selected  features,  XBT 
depths  calculated  with  coefficients  resulting  from  modified  overall 
regression  (ZXbt  =  6.798/  -  2.38  X  10~¥  -  4.02).  Different  symbols 
denote  different  LAS’s. 


What  is  quite  apparent  is  that  a  great  deal  of  variability 
remains  from  probe  to  probe,  although  the  systematic 
skewing  of  depth  errors  (Fig.  4)  has  been  removed. 

c.  Individual  profile  regressions  and  corrections 

Toward  a  better  understanding  of  the  variability  of 
depth  errors,  we  conducted  a  modified  regression  anal¬ 
ysis  (with  7?  =  2.06  X  1(T4  m_1),  producing  a  set  of 
coefficients  for  each  profile.  Subsequently,  the  depths 
in  each  profile  were  corrected  using  its  respective  set 
of  coefficients,  and  depth  differences  were  again  plotted 
( Fig.  7 );  a  significant  reduction  in  the  errors  is  apparent, 
but  a  slight  trend  remains  in  the  deeper  layers. 

Analogous  to  Figs.  2a, b,  the  same  set  of  simultaneous 
XBTs,  with  depths  corrected  with  individually  deter¬ 
mined  coefficients,  is  plotted  in  Figs.  8a,b.  This  is  typ¬ 
ical  of  all  individually  corrected  XBTs  and  illustrates 
the  empirical  validity  of  the  form  of  the  FRE  used.  We 
conclude  here  that  most  of  the  variability  from  probe 
to  probe  is  accounted  for  by  variations  in  the  linear 
and  constant  terms,  and  not  in  r\\  that  one  equation  of 
the  form  ( 2 1 )  is  adequate  to  represent  the  fall  rate  of 
the  XBT  for  depths  >  10  m;  and  that  more  complex 
equations  involving  higher-order  terms  are  probably 
unnecessary. 

It  is  useful  to  examine  more  closely  the  variability 
of  the  coefficients.  Questions  arise  as  to  whether  most 
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Fig  .  7.  Depth  error  (Zero  -  ZXBt  )  for  selected  features  and  XBT 
d  pths  calculated  using  individually  determined  coefficients.  Different 
symbols  denote  different  LAS’s. 


of  the  variability  is  temporal,  the  result  of  differences 
in  LAS  problems,  or  just  due  to  differences  in  the 
probes  themselves.  The  coefficients  ( b  and  c)  found 
using  i]  =  2.06  X  10 ~3  m-1  are  plotted  as  functions  of 
time  in  Figs.  9  and  10.  The  different  symbols  in  the 
figures  denote  which  LAS  is  represented.  There  is  no 
clear  pattern  to  the  variations  in  b  and  c.  They  are  not 
correlated,  and  there  appears  to  be  as  much  variation 
between  values  for  fixed  times  as  there  is  for  a  given 
LAS  as  a  function  of  time.  For  this  dataset,  much  of 
the  scatter  in  c  can  be  attributed  to  a  LAS  malfunction 
(Sippican,  personal  communication)  where  the  first 
few  scans  of  data  (the  exact  number  is  unknown  and 
appears  to  be  random)  were  lost  at  the  beginning  of 
each  probe  launch.  Means  and  standard  deviations  of 
b  and  c  quantities  appear  in  Figs.  1 1  and  12.  The  right 
ordinate  axis  of  Fig.  11  is  in  units  of  depth  and  rep¬ 
resents  the  contribution  of  the  linear  term  to  the  FRE 
at  /  =  1 22  s,  the  time  in  which  the  probe  reaches  about 
760  m  using  the  original  Sippican  coefficients.  Also 
appearing  is  the  rms  spread,  based  on  the  maximum 
difference  of  b  (and  of  c)  within  each  simultaneous 
XBT  group.  There  is  no  LAS  dependence  apparent  for 
c,  so  the  acquisition  malfunction  appears  to  have  been 
ubiquitous.  It  should  be  noted  here  that  XCTD  probes 
launched  from  these  same  systems  on  the  same  cruise 
exhibited  standard  deviations  in  c  of  only  about  1 .5  m 
( Hallock  and  Teague  1 990 ) .  For  the  XCTDs,  the  LAS’s 


Temperature(deg  C) 


Fig.  8.  Typical  set  of  simultaneous  XBT  profiles  (solid  curves) 
with  depths  calculated  with  individually  determined  coefficients. 
Time-interpolated  CTD  temperature  (dashed)  is  also  plotted:  (a) 
upper  140  m,  (b)  400-700  m. 
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FIG.  9.  Individually  determined  coefficient  b.  Different  symbols  denote  different  LAS’s. 


were  in  a  different  mode  and  data  were  not  lost  (Sip- 
pican,  personal  communication).  Standard  deviations 
and  rms  spread  are  virtually  the  same  for  b  (and  for 
c).  The  mean  value  of  b  for  LAS  3  is  barely  significantly 
higher  than  for  the  other  LAS’s.  Since  there  were  only 
18  XBT  profiles  acquired  with  LAS  3  (versus  about 
33  for  the  others),  we  suspected  a  bias  due  to  the  miss¬ 
ing  data.  We  recalculated  the  means  of  data  acquired 
with  LAS  1 ,  LAS  2,  and  LAS  4  including  only  profiles 
corresponding  to  those  from  LAS  3;  there  were  virtually 
no  differences  from  the  values  appearing  in  Fig.  11. 
The  reason  for  the  anomalous  mean  b  for  LAS  3  is  not 
clear. 

For  comparison,  b  from  several  other  studies,  as  well 
as  the  Sippican  value,  also  appears  in  Fig.  1 1.  The  value 
found  by  Singer  ( 1990)  is  within  the  standard  deviation 
envelopes  for  three  of  the  four  records  of  the  present 
study.  That  of  Hanawa  and  Yoritaka  ( 1987)  is  barely 
outside  the  envelopes.  This  suggests  that  hypothesized 
regional  differences  in  fall  rates  may  not  be  statistically 
significant. 


A  portion  of  the  variability  in  c  (and  perhaps  part 
of  the  mean  c  as  well )  is  the  result  of  a  data-acquisition 
problem.  The  remainder  may  well  be  due  to  acceler¬ 
ation  of  the  probes  within  the  first  10  m;  differences 
in  probe  entry  angle  or  other  mechanical  perturbations 
may  contribute  to  the  observed  scatter.  Unfortunately, 
we  cannot  separate  the  LAS  problem  from  other  effects 
with  this  dataset,  so  a  definitive  determination  of  a  real 
offset  for  use  in  Eq.  ( 15)  awaits  a  future  investigation. 
The  primary  cause  of  probe-to-probe  variability  is  the 
coefficient  b  (or  WT).  While  overall  or  average  values 
of  b  may  ary  due  to  environmental  conditions  (e.g., 
the  density  profile),  probe-to-probe  variability  seen  in 
the  present  study  is  most  likely  the  result  of  differing 
probe  weight,  volume  and/or  effective  drag.  G84  con¬ 
cluded  (based  on  actual  measurements)  that  probe 
weight  and  volume  probably  cannot  explain  the  vari¬ 
ability,  implying  that  differences  in  drag  are  the  cause. 
These  might  include  slight  differences  in  roughness  of 
the  surfaces  of  the  probes,  affecting  the  separation  point 
for  turbulent  layers. 
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Fig  .  10.  Individually  determined  coefficient  c.  Different  symbols  denote  different  LAS’s. 
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Fig.  1 1.  Means  (dots)  and  standard  deviations  (error  bars)  of 
individually  determined  regression  coefficients  b  for  the  four  LAS's 
(labeled  1.  2.  3.  4):  rms  of  maximum  spread  of  b  in  simultaneous 
X  BT  groups  ( labeled  sp ) .  Also  shown  are  values  of  h  found  by  other 
investigators:  WS — Wright  and  Szabados  ( 1987 ),  Sg — Singer  ( 1990). 
HY — Hanawa  and  Yoritaka  ( 1989),  Sp — original  Sippican  value. 


d.  Implications  for  the  entire  temperature  profile 

The  preceding  analyses  and  dicussion  were  based  on 
selected  points  of  the  temperature  profile.  These  were 
chosen  to  avoid  contamination  of  results  by  errors  in 
measured  temperature.  The  depth  interval  between 
about  1 50  m  and  400  m  and  that  deeper  than  about 
630  m  were  not  directly  addressed.  Hence,  it  was  useful 
to  assess  the  effectiveness  of  the  new  FRE  ( second  set 
of  coefficients  in  Table  2)  in  accurately  representing 


Fig.  12.  Means  (dots)  and  standard  deviations  (error  bars)  of 
individually  determined  regression  coefficients  c  for  the  four  LASs 
(labeled  I.  2,  3,  4);  rms  of  maximum  spread  of  c  in  simultaneous 
XBT  groups  (labeled  sp). 
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the  temperature  profile  over  its  entire  depth  range  ( 10- 
760  m).  To  do  this  we  formed  an  error  profile 

5T(Z)  =  Tctd(Z)  -  Txbt(Z)  (23) 

for  each  XBT  profile.  Here  7'CTD  is  CTD  temperature 
data  interpolated  to  the  time  of  the  corresponding  XBT 
sample  (explained  earlier).  We  then  calculated  mean 
and  rms  profiles  by  averaging  over  all  118  XBTs.  For 
reference  purposes  we  first  examine  the  statistics  of  bT 
calculated  using  the  uncorrected  XBT  profiles.  Figure 
1 3  (analogous  to  Fig.  4  for  the  depths  of  features )  shows 
a  mean-error  profile  that  starts  near  -0.2°  near  the 
surface,  increases  to  a  maximum  near  0.4°,  and  de¬ 
creases  again  to  about  0.1 2°C.  Superimposed  on  this 
overall  trend  are  excursions  of  the  order  of  0.2°C  over 
vertical  scales  of  about  20  m.  The  rms  error  profile 
reflects  primarily  the  mean.  By  recalculating  XBT 
depths  using  the  new  FRE.  we  have  removed  most  of 
the  large-scale  trend  in  bT  (Fig.  13;  analogous  to  Fig. 
6)  and  have  reduced  the  rms  error  to  a  considerable 
degree.  The  remaining  error  profile  has  two  aspects  of 
note;  the  same  ( but  somewhat  reduced )  vertical  struc¬ 
ture  near  the  surface  and  in  the  depth  aperture  con¬ 
taining  the  step  features,  and  a  nearly  depth-indepen¬ 
dent  negative  offset  of  the  mean  error.  We  suggest  that 


Temperature  Error  (deg  C) 
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Fig.  13.  Temperature  error  (  7V-rD  -  TXbt  )  mean  (solid) and  stan¬ 
dard  deviation  (dashed)  profiles  for  XBT  depths  calculated  with  the 
Sippican  FRE. 
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the  approximate  20-m  vertical  structure  in  the  error 
profile  is  the  result  of  probe-to-probe  variability  in  fall- 
rate  parameters  as  discussed  in  the  preceding  subsec¬ 
tion.  This  hypothesis  is  addressed  in  the  following.  The 
nearly  constant  offset  error  cannot  be  the  result  of  a 
residual  error  in  the  FRE  coefficients  since  the  tem¬ 
perature  profile  is  not  a  linear  function  of  depth.  Ex¬ 
amination  of  individual  XBT  error  profiles  (not 
shown )  showed  them  to  be  primarily  depth-indepen- 
dent  offsets.  The  mean  offset  (Fig.  1 4 )  is  about  -0. 1  °C. 
This  translates  to  an  isotherm  depth  error,  according 
to 

5Z  =  (dT/dZy'6T.  (24) 

Near  the  bottom  of  the  profile,  dT/dZ  «=  1 0-2 °C  m~ 1 , 
yielding  6Z  »  10  m.  Individual  offsets  can  be  nearly 
twice  this  value. 

'  To  eliminate  any  ambiguity  as  to  the  source  of  the 
temperature  offset  error,  we  used  an  approach  analo¬ 
gous  to  that  used  to  calculate  the  FRE  coefficients. 
Here,  however,  we  identified  eight  points  in  the  CTD 
temperature  profiles  that  were  relatively  isothermal  in 
depth  and  time.  The  characteristics  of  these  features 
are  summarized  in  Table  3.  Six  of  the  points  were  in 
the  virtually  isothermal  layers  below  the  steps  used  in 


Temperature  Error  (deg  C) 

-0.6  -0.4  -0.2  0.0  0.2  0.4  0.6 


Fig.  14.  Temperature  error (Tctd-  7"xbt)  mean  (solid)  and  stan¬ 
dard  deviation  ( dashed )  profiles  for  XBT  depths  calculated  with  the 
second  FRE  of  Table  2  [Eq.  (25)]. 


Table  3.  Selected  temperature  calibration  points. 


Z~td  (m) 

rcTD(°C) 

^co(°C) 

8 

27.64 

0.02 

83 

25.83 

0.01 

446 

11.66 

0.00 

507 

10.32 

0.01 

530 

9.69 

0.00 

551 

9.22 

0.01 

572 

8.71 

0.01 

595 

8.35 

0.02 

the  previous  (FRE)  analyses.  The  remaining  two  points 
were  in  the  upper  100  m.  The  standard  deviations  for 
all  selected  points  are  less  than  0.03°C.  The  same  fea¬ 
tures  were  identified  in  the  XBT  profiles,  and  profiles 
of  AT  [defined  as  in  Eq.  (23)  but  for  the  selected  points] 
were  derived.  The  vertical  gradients  of  XBT  temper¬ 
ature  for  the  selected  points  were  at  most  0.03°C  m- 1 . 
We  reemphasize  here  that  AT  found  in  this  way  is  vir¬ 
tually  independent  of  FRE  errors.  We  then  calculated 
an  average  (over  depth)  AT  for  each  XBT  and  cor¬ 
rected  each  actual  XBT  profile  with  this  average.  Figure 
1 5  depicts  the  error  profile  statistics  following  temper- 


Temperature  Error  (deg  C) 
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Fig.  15.  Temperature  error  (Tctd-  Txrr)  mean  (solid)  and  stan¬ 
dard  deviation  (dashed)  profiles  for  XBT  depths  calculated  with  the 
second  FRE  of  Table  2  [Eq.  (25)]  and  with  corrected  XBT  temper¬ 
atures. 
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ature  correction  (using  the  XBTs  with  depths  calcu¬ 
lated  with  the  new  FRE ) .  Most  of  the  offset  error  has 
been  removed,  and  the  rms  error  has  been  reduced 
somewhat  below  about  550  m. 

Finally,  to  test  the  assertion  made  earlier  than  most 
of  the  remaining  structure  and  variability  in  the  error 
profile  are  the  result  of  probe-to-probe  variations  in 
fall-rate  parameters  we  further  depth  corrected  each 
XBT  with  its  respective  individually  determined  FRE 
coefficients  (as  in  the  previous  subsection).  The  re¬ 
sulting  error  profiles  appear  in  Fig.  16  (analogous  to 
Fig.  7).  Indeed,  most  of  the  structure  and  variability 
are  gone.  The  largest  remaining  errors  occur  between 
100  and  250  m  and  are  generally  less  than  0. 1  °C. 

In  summary,  the  proposed  new  FRE  coefficients  ef¬ 
fect  a  significant  improvement  in  temperature  mea¬ 
surement  capability  over  the  full  depth  range  of  the 
T-7  XBT.  Most  of  the  remaining  errors  are  due  to  in¬ 
dividual  physical  probe  differences  and  depth-inde¬ 
pendent  temperature  bias.  A  bulk,  quantitative  mea¬ 
sure  of  the  improvement  is  given  in  Table  4,  which 
shows  5T statistics  computed  over  all  XBTs  and  depth. 
The  largest  improvement  is  the  decrease  in  the  mean 
error,  while  the  rms  error  was  reduced  by  about  40%. 


Temperature  Error  (deg  C) 


Fig  .  1 6.  Temperature  error  (  Tctd  -  7x bt  )  mean  ( solid )  and  stan¬ 
dard  deviation  (dashed)  profiles  for  XBT  depths  calculated  using 
individually  determined  FREs  (section  4c)  and  with  corrected  XBT 
temperatures. 


Table  4.  Overall  temperature  errors. 

FRE 

5T(°  C) 

rms  (ST)  (°C) 

Sippican 

0.19 

0.32 

Proposed 

0.07 

0.19 

5.  Conclusions  and  recommendations 

With  the  analysis  of  concurrent  CTD  and  XBT  tem¬ 
perature  profile  data  acquired  in  the  Atlantic  Ocean 
near  Barbados  we  have  shown  (as  have  others)  that 
the  presently  used  T-7  fall-rate  equation  is  inadequate. 
An  improved  equation  is  proposed  for  depths  greater 
than  10  m: 

Z  =  6.798/ -  0.00238?2  -  4.01.  (25) 

This  new  equation  reduces  overall  rms  XBT-measured 
temperature  error  by  about  30%.  The  rms  residual 
depth  error  in  the  determination  of  the  coefficients  of 
Eq.  ( 25 )  is  about  5.2  m.  Most  of  this  error  is  the  result 
of  probe-to-probe  differences  in  terminal  velocities, 
which,  in  turn,  are  probably  due  to  variations  in  probe- 
drag  characteristics.  The  probe-to-probe  variability  di¬ 
minishes  the  relative  importance  of  differences  in  fall 
rates  found  in  different  studies,  and  is  a  caveat  to  in¬ 
vestigators  who  use  XBTs  to  examine  finescale  features. 
The  remedy  for  the  variability  will  first  require  deter¬ 
mination  of  its  cause.  If  it  is  the  result  of  differences 
in  effective  drag,  some  way  must  be  found  to  make  the 
probes  more  uniform  in  this  respect,  perhaps  involving 
mechanical  design  modifications. 

We  further  propose  a  more  complete  equation  that 
includes  an  initial  acceleration  in  the  upper  10  m, 

Z  =  Zi(t,  IT0)  -  0.00238/2,  (26) 

where  Zt  [Eq.  ( 14)]  accounts  for  the  near-surface  ac¬ 
celeration  and  approaches  the  linear  and  constant  terms 
of  Eq.  (25)  as  depth  approaches  10  m.  Parameter  W0 
is  the  initial  downward  velocity  at  the  surface.  The 
present  dataset  is  inadequate  to  determine  accurately 
the  parameters  of  Z\ ,  but  a  value  of  W0  =  0  would 
yield  empirically  the  offset  of  about  -4  m  found  in  the 
analysis.  Further  study  of  the  near-surface  transient 
problem  is  indicated.  In  particular,  an  experiment  in 
water  at  least  1 0  m  deep  will  be  required  where  probes 
can  be  tracked  to  within  a  few  centimeters  and  timed 
to  within  a  few  milliseconds.  During  this  experiment, 
several  series  of  probe  drops  should  be  made,  varying 
parameters  such  as  height  above  the  surface,  angle  from 
the  vertical,  and  amount  of  vertical  rotation  (i.e., 
“wobble”).  Spinup  (about  probe  axis)  is  also  likely  to 
be  a  factor  in  the  overall  acceleration  and  should  be 
addressed. 
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